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ABSTRACT 


Finite-amplitude standing waves in air at ambient con- 
Olurrons COntained in a rigid-walled cylindrical ‘ele with 
a large length-to-diameter ratio were experimentally inves- 
tigated. The pressure waveform at the end of the tube was 
digitized and Fourier analyzed on an IBM 360 digital 
computer. Amplitudes and phases were obtained for all 
harmonics with amplitudes greater than 1% of the fundamental 
mest benign) parameters from 0.25 to 1.90 and for frequency 
parameters from -0.8 to 2.0. The strength and geen 
parameters are defined as MbO and 2hE/Ot, respectively, 
Miers tne reen numoer Of fundamental, b the nonlin- 
Habmevyeparameter, © the quality factor of the resonator, 
Af the frequency away from fundamental resonance, and 


feeeeene band width at the half-power points. When these 


= 
results are compared to the theoretical model of Coppens 
and Sanders, it is seen that while the theory accurately 
predicts the magnitude and shape of the harmonic content, 
1t consistently underestimates the AP @gREEINeN at which each 


harmonic peaks. In addition, the theory fails to predict, 


except qualitatively, the phase angles. 
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Pree LE RODUCTION 


Coppens and Sanders [1,2] have developed a one- 
dimensional, nonlinear wave equation with dissipative terms 
corresponding to those encountered in a rigid-walled closed 
tube with a large length-to-diameter ratio. An attenuation 


@enstant of the form 


is assumed, where a. Meee ne absorption Bee ioiene associa- 
ted with the frequency nw which characterizes the wall 
losses, Ris the radius of the tube and L is the length of 
the tube. With this additional loss mechanism the wave 


equation takes the form 














. ae at + D_| UL = ie o§ 
il yx2 eee 542 Te 1h Oxo \ ox} © 
where 
eel) 2 (C,/C41)/2, 
a 1 52) eee 
SMe cio oe * 3x7} 
% 
$1 = & lao fis Al 
yY = oy the ratio of specific heats, 
v = kinematic viscosity, 
Pr = Prandtl number, 
W = angular frequency of the fundamental, 


Den harmonic number, 





Ces the thermodynamic speed of sound, 
a, = absorption coefficient for the lowest frequency 
and 
. Mes 
U = particle velocity = ye se #1) 
n 


ae 


To arrive at the above equation it was necessary to make the 

simplifying assumptions that the wall losses were much 

greater than the bulk losses and that the Mach number 

M= P./p C - 
Lo) Me, 


is much less than unity. Py is the peak value of the funda- 
mental of the pressure at the rigid end of the tube. 
A Fourier series solution of the form 


ee (sents) 
fe ies ee 


co 
vw" \ 
ove, &(n) Cos int (1- 7,)| Sin (nwtt+on) 


1 ie 
1s assumed with the tube excited with a pressure distribu- 


tion given by 
Z ; 
Py (x,t) = aes Gost Gx i) Sinwt . 


(tei iomool MGlon was Substituted directly into the non- 
linear wave equation, the following set of algebraic 
equations resulted: 


H [cose Noun (sins )c,| 
n nN n n n 


PeutsumnG.+¢...5.)- 2 (Ss ) 3 


ee ee ae 


and 





H (cose, }C,, i (Sine, )s_ | 


1 = IL 
= {= £ (Cc . - S Se 
2% 2 5=1 Oe ae N=) ° 3 
- 2 (Con n+ SM}, 
j=1 a) a) 
where 
Sy eRe oni O. , Ch = RL Cos dnt 
d = 2/SP, 
-; 
-] I oa = 2Aw/w 64 
a = Tan = _ 
ie 
n tna,/o, 
aol 
= dhishel [A/B], 
%5 
H = Tees ; 
n 
a, = infinitesimal amplitude absorption coefficient 
in free space, 
wo. = angular frequency at resonance, 
SP = strength parameter 
_ _ 2 
SP = Mbk/2 i= P bkK/ (2p .C. a,) 
k = wave number 
and 
Aw = WwW - w_. 
Ie 


These equations were then programmed for computer solution. 
For given values of Aw and SP, the computer calculated 
values of Pi/Py and Oo. At present, solutions have been 


found for strength parameters up to 0.75. 





Previous investigations in this area have heen carried 
out by Sanders [1,2] and Beech [3]. Both investigators 
examined the harmonic content of the pressure waveform at 
the rigid end of the tube near resonance conditions. 
Sanders employed a wave analyzer to obtain the frequency 
response for the first four overtones. Good agreement with 
theory was found but no phase information was obtained. 
Beech used photographic and graphical techniques to provide 
numerical analysis of pressure waveforms near resonance. 
Phase MEGrHAr Ton, however, was obtained only for cases at 
resonance. Again good agreement with theory was obtained 
for cases approaching the onset of shock. 

The present investigation examines both the harmonic 
content and phase relationships for strength parameters as 
high as 0.75. Again numerical methods were employed. 
Whereas the graphical techniques of Beech were used to 
generate 64 data equally spaced throughout the pressure 
wave form, the technique in this investigation employed 
analog-to-digital conversion by computer to provide better 
than 1000 data to represent the pressure waveform. With 
the pressure waveform digitized to such Fe Gee it was 
possible to obtain amplitude and phase information for all 
harmonics with amplitudes greater than one percent of the 


fumeaamnental . 





Pie oe IMENTAL CONSIDERATIONS 


fae APPARATUS 

A block diagram of the experimental system used in this 
investigation is shown in Figure 1. Finite-ampiitude 
standing waves were generated in air at ambient conditions 
contained within a tube six-feet long with an inner diameter 
eee. 250 an. and a wall thickness of 1.125 in. The air was 
excited at one end by a piston which was driven by an M-B 
Electronics Model EA 1500 exciter capable of a maximum 
no-load acceleration of 124g. The FA 1500 was Been by two 
M-B electronics Model 2120MB power amplifiers operating in 
Pavtallei. The Waxemum attainable acceleration with the 
Beugqhet Of the piston as load was about 50g. The frequency 
and power level were controlled by a General Radio Type 
1161-AGC coherent decade frequency synthesizer. With this 
synthesizer it was possible to select frequency to better 
than 0.001 cycles/sec. (Monitoring the frequency with a 
Hewlett-Packard Model 521C frequency counter showed the 
frequency stability to be better than one paras alia, ING) aig ealaveyer,) 
The acceleration of the piston was determined by measuring 
the output an Endeveco Model 2215 accelerometer implanted 
within the piston. When the piston was withdrawn from the 
abe and driven in free air at the approximate experimental 
frequency and acceleration, the harmonic distortion was 


measured to be less than 0.3 percent. The other end of tne 





tube was capped with a thick plate which was securely bolted 
to the tube. The sound pressure level at the rigid end of 
the tube was detected with a 1/4=-in. diameter Bruel and 
Kjaer Type 4136 condenser microphone mounted in the cap so 
that its diaphragm was flush with the end ofthe tube. The 
iverspnone response was flat to within + 0.5 dB from 80 Hz 
to 20 kHz. The outputs of the microphone and accelerometer 
were monitored with Hewlett-Packard 400D vacuum-tuhbe volt- 
meters. The microphone output also went to a Hewlett-Packard 
Model 302A wave analyzer so that the first two overtones 
could be monitored. The microphone output was recorded by 
a PI 6200 General Purpose Portable TS crementacion Tape 
Recorder. Recordings were made in the FM mode at a tape 
Breed or 37.5 in/s so that the froaquency response of the 
mapas eecOrdc,r was +t 1 dqdB from DC to 10 kHz. Although the 
Signal was recorded at 37.5 in/s it was necessary when 
Seer zing the signal to play back at a speed of 3.75 in/s. 
To insure that this change in tape speeds did not alter 

the recorded frequencies, a test Saba was recorded at the 
higher speed and then played back at the slow speed. The 
two signals were observed on a dual trace seal legen 
where it was observed that the frequency of the test signal 
was 10.0 times that of the signal played back at 3.75 in. 


per second with no discernable drift between the two signals. 


B. MICROPHONE CALIBRATION 
Investigation of the standing waves at various strength 


es, 


parameters requires an absolute calibration of the microphone 


10 





system. The microphone system is defined here as the micro- 
phone, preamplifier and the equipment loaded onto the output 


of the preamplifier. The system sensitivity is defined as 


S, = V,/? 


where V_, is the peak voltage output of the preamplifier and 
P ts the peak pressure at the microphone. Two methods were 
employed to determine S,. The first and simplest used a 
Bruel and Kjaer Model 4220 pistonphone which produces a 
known sound pressure level at 250 Hz. The sensitivity was 
found to be Sm wae eet Ol. 2) AB eed VOLEZ tbat, where the 
uncertainty ie that inherent in the pistonphone specifications. 
The second method makes use of theoretical predictions by 
Coppens for small strength parameters as verified by Beech 
[3]. Coppens shows that at resonance the ratio of the peak 
pressure of the second harmonic to the square of the peak 


pressure of the fundamental 1s given by 


P,/Ps = oi wB)/ (89 Coa) . 


When the definition of sensitivity is substituted into the 
above equation the result is 

S. = (Bwv7)/ (47 0 ea ) 

m il OO 


Sin was then found by exciting the tube at resonance and 
measuring the fundamental and second-harmonic voltages. By 
this method = Steve Ons) AB xe ) volt/ubar where the 
uncertainty is due largely to the error in determining the 


Sereomuatlon constant. 


deh 





SeeeMicROPHONE OUTPUT AND STRENGTH PARAMFTER 
In this investigation, where the pressure waveform was 
to be analyzed at various strength parameters, it was 


necessary to find a relationship between the strength 








parameter and the observable quantities. Using the 
equations 
eee log|—-=-| + 72.0 
me 77S es 
Pp ly im 
ee = spr, - 197.0; 
2 ES C 2 
Povo Povo 
and 
M = SP(20/K) 
b 


Geemme@an Show that the ris value of the fundamental cf the 
microphone output (v) is related to strength parameter and 


the attenuation constant by 
WineeL SOrnWeCL 2 oP 
The constant 136.0 is relatively insensitive to small 


temperature variations about room temperature. 


D. ATTENUATION CONSTANT 
The attenuation constant for infinitesimal amplitudes 


has been shown by Beech to be 


AMIE 








le 
Cc 

O 
where Af is the frequency difference between the half power 


boilmtes On the resonance curve (taken for constant accelera- 


tion amplitude of the piston). 
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Peer REOURNCY PARAMETER 


The frequency parameter is defined by 


FP = 2hw/ 554 
where 

a= Wb = W 
and 

4 = 20,C /W- 


Mmreserrequcncy parameter 1s related to the experimental 


Frequency increment (Af...) by 


Af ‘= Af + FP/2, 
el. 


where Af is the frequency difference between the half power 
pemiges OL the infijmtesimal resonance curve as used in 
calculating Oy in this investigation Af was found to be 


eset 0.015 so that 


3 


Af = Oe On iP 
ex 


me) 





viet COMPUPAR ANALYSIS 


fe OAMPLE RATE 

Many computer programs exist which, when supplied with 
data points defining a periodic waveform, will give the 
Fourier coefficients. The first consideration, generally, 
is the number of data required per cycle. Theoretically 
this number should be at least twice the number of harmonics 
desired. However, the question of accuracy of the computer 
predictions of harmonic amplitudes arises immediately. As 
might be expected, the accuracy of the harmonics Beearmined 
tails off faster as fewer samples per cycle are used. This 
Pee ddsrtoeetne problem of first datcrmining how many 
Significant harmonics are desired. For finite-amplitude 
standing waves of large strength parameter the waveform 
Beeeetatively approaches a function whose harmonics content 
Falls off inversely with harmonic number. If such a wave- 
form is passed through amplifiers Bee torage devices before 
being analyzed, harmonics below a given amplitude will he 
loc tmanetae simnierent moilse. The device i Sane investiga- 
tion which was most seriously restrictive in this respect 
was the PI 6200 FM tape recorder. It was found that the 
dynamic range of this instrument was about 42 db which 
restricted the investigation to harmonics with amplitudes 
greater than about one percent of the fundamental. Fora 
Dowteclewaamp function this restriction occurs at about the 


Btn Marmon ic. 
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Determination of the number of samples per cycle to be 
used was based on the arbitrary criterion that the error 
in the 30th harmonic be less than 0.1 percent for an ideal 
ramp function, Applying several computer algorithms for 
determining Fourier coefficients to an ideal ramp function 
one finds that 1000 samples per cycle are required for the 


above accuracy. 


B. THE DIGITAL-COMPUTER PROGRAM 

The primary function of this program was to generate 
Seamer zcd HOurier coefficients for a periodic function 
Given digitized data on a seven-track tape. Secondary to 
this, the program changed the form of the coefficients to 
A(n) Sin(Nwt+  o(n) and computed A(n)/A(1) for each harmonic. 
Optional functions provided by the program were the averag- 
ing of several periods of data to reduce the affects of 
Gandom noise. An optional plot of the input data was avail- 
able along with points reconstructed from the calculated 
Peou@lctecoecrttteients. This plot provided a visual check on 
the input data and the computed coefficients. 

The digitized data provided by the SDS.9300 were written 
On seven-track tape. HONSWeE the binary representation of 
numbers on the seven-track tape differs from those on the 
standard nine-track tape used by the IBM 360 system. To 
accomodate the seven-track tape, a conversion program 
(subroutine FORM) was used to convert the data. The con- 
Voeee-eecaca weremnot floating point but integers and required 


further conversion with the use of a scaling factor. 
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Prime consideration was given to the types of algorithms 
commonly used in generating Fourier coefficients. Two of 
the more widely used methods are the Seo lley tulysy [5] Fast 
Fourier transform and the recursive technique of Goertzel 
[4]. Both methods were tested for accuracy in predicting 
the coefficients of known waveforms. The two methods were 
found to be equal. One advantage of the Cooey-Tukey method 
is its minimal use of computer time. Its computation time 
is about an order of magnitude less than that required by 
the Goertzel method. However, when this time saving is 
compared with the overall required computer time for analysis 
the advantage diminishes. The final sone eration eValel 4) ale 
the end the most important, was the flexibility in using 
Mom men CcoCley-Turay mMEEnod required an intecer sower-of- 
two number of points. To he able to satisfy this condition 
along with the requirement for 1000 or more samples per cycle 
complicated the digitizing process. The Goertzel method 
required only an odd number of samples and for this reason 


was used in the computer program. 


Ce LHF ANALOG TO DIGITAL CONVERSION (A-D) . 

The conversion of the eee data was accomplished with 
Pisco wOtucd ivemra computer consisting of the SNS 9300 
high-speed digital computer and the Ci 5000 general-purpose 
analog computer. |Use of Ehe system for A-D conversion 
requires a FORTRAN program for the SDS 9300 and properly 
Senstructed analog and logical patch boards for the Ci 50CO. 


In addition, several peripheral equipments are required. 
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These are the card reader, magnetic tape transport, line 
printer and, if desired, a pen recorder. 

Aside from the question of sampling, which was deter- 
mined earlier, the important question in the A-D conversion 
is accuracy. The clock in the computer is accurate to one 
Dart in 100 million, well within any accuracy requirements. 
The A-D converter is designed to accept signals between 
fee volts and divide this interval into 16,384 parts 
Beevtaing a resolution of 12.2 mV. For a properly amplified 
analog signal, the dynamic range of the A-D converter is 
much greater than the PI 6200 tape recorder. 

The analog patch board (Figure 3), as used one Pee CO 
version, served to amplify the analog signal prior to con- 
Twemeremo and Gasteibucea the signal to the converter and 
oscilliscope. Additionally, a patch was provided linking 
the Digital-to-Analog converter with a pen recorder. This 
additional function is explained in detail in the pro- 
cedures section. 

The logic patch board (Figure 4), had three functions. 
The first provided a sampling frequency (10,000 Hz) to the 
A-D converter. This frequency represents oe approximate 
maximum sampling rate allowed. Since the approximate 
yesonant frequency of the tube was 94 Hz it was necessary to 
play back signals from the PI 6200 at one tenth the recorded 
speed Pomacauure the desirea number of samples per cycle. 
The second function provided a triggering voltage through 


pie sternal switch (DS-1) which was used to start the data 


i) 





conversion process. Associated with this function is an 
externally controlled delay (DF00) which insured that the 
triggering voltage remained on long enough for one set of 
data (one record) to be written on the seven-track tape. 

The last function provided another triggering voltage 
through an external switch (DS-2). When the last record was 
to be written on tape switch DS-2 was turned on while 


Switch DS-2 was triggered. 


PeoeeorvolhM TEST 

Prior to analyzing experimental data, it was necessary 
to test the entire system process of A-D conversion and 
computer analysis. The results of this test for a ramp 


Pomecton 1S shown in Figure 5. 
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IV. PROCFDURFS 


Poeeee RERUN PROCEDURES 

Beech pointed out that the resonant frequency of the 
tube changed by 0.16 Hz for each degree centigrade of 
temperature change. To minimize this affect, all data 
gathering runs were made during the evening hours when 
the room temperature remained relatively constant for 
several hours. Typical resonant frequency variations during 
these hours were of the order of 0.01 Hz per hour. Several 
hours prior to taking data the equipment was Eee on 
and allowed to stabilize. During this warmup period, the 
Peeper Womcape Decorcer waS Calibrated tc insuxe that the 
FM system was properly aligned and to insure that the 
Signals being recorded were not overdriving the tape recor- 
der preamplifier and thereby introducing unwanted signal 
distortion. Next, the alignment of the piston in the free 
end of the tube was checked by observing the output of the 
accelerometer on an oscilloscope and measuring, with a 
wave analyzer, the second- and ind ed herenerane Content wt 
necessary, the harmonic content could be reduced by carefully 
adjusting the alignment of the piston so that the distortion 
from any harmonic was less than 2 percent. To Minimize 
leakage of sound from the tube, the piston was fitted with 
oO itngeana lubricated. The harmonic distortion intro- 


duced by the piston was sensitive to the type of lubricant 
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used. It was found that a light turbine oil minimized the 
Gistortion while providing an adequate seal. 

Just prior to taking each set of data, it was necessary 
to determine the resonant frequency and the attenuation 
constant. For this measurement the piston was excited at 
a constant value of acceleration low enough so that the 
strength parameter at resonance was lessthan 0.03. From 
the frequency difference between the half power points, the 
attenuation constant was determined. The resonance frequency 


was taken to lie half way between the half power points. 


EeeeeATA TAKING PROCEDURE 

Corresponding to the desired strength parameter was 
a given value of the pressure fundamental as shown earlier. 
With the piston excitation set to the proper level to keep 
the pressure fundamental constant, data were recorded at 
O.1 Hz intervals for twenty frequencies from 0.6 Hz below 
resonance to 1.3 Hz above. At each frequency, the signals 
from the microphone and accelerometer were recorded on the 
Peo O0matehi ga speed (37.5 in/s) for about 10 seconds 
(30 £ft.), and the second-harmonic amplitude of the pressure 
was noted to be later mesated with the computer analysis. 
Since the taking of an entire set of data at one strength 
parameter required from 30 minutes to an hour, at the end 
of the run it was necessary to check the resonant frequency 
and attenuation constant to insure that the variation in 
resonant frequency was substantially less than one frequency 
iieervaleG.ileHz) and that the change in attenuation constant 
was less than one percent. 


Z0 








The theoretical predictions (8) are compared with the 
experimental results in Figures 6 through 13 for each value 
of strength parameter, plots are shown of the harmonic 
amplitude vs. frequency parameter and of the phase difference 
between consecutive harmonics vs. inverse harmonic number. 
Figure 6 shows the actual data obtained for three different 
experimental runs. The degree of scatter is a demonstra- 
tion of the excellent reproducibility of the results. Both 
experimental and theoretical results are limited to harmonics 
with amplitude greater than 1% of the fundamental. 

In comparing the experimental values of the harmonic. 
amplitudes with the theoretical predictions for a strength 
parameter of 0.25 it is seen that the agreement is fairly 
good but that there does seem tc he a slight systemat c 
difference. The corresponding curves for higher strength 
parameters display this difference to a more marked degree. 
It is seen that the experimental results peak at about the 


predicted values of Piny/? ¢ and that the shape of the 


1) 
curve for each harmonic is about as predicted. However, 
the experimental Bees Weak at a value of frequency 
parameter significantly greater than predicted and this 
difference increases with harmonic number. It is interest- 
ing to note that at the frequency of maximum distortion of 
the pressure waveform (at approximately FP = 0.5) theory 
Accurately predicts the harmonic amplitudes. 


The comparison between the experimental and theoretical 


phase angles displays an even greater discrepancy. While 





values of o, - eT do seem to be linearly related to l/n, 
the values for the second harmonic are greatly different 
and the slopes of the experimental curves are significantly 
greater than predicted. 

In conclusion, it appears that while the present theo- 
retical model is able to predict the magnitude and shape of 
mpenmnarmonic content of the distorted waveform, it consis- 
tently underestimates the frequency at which each harmonic 
peaks. This underestimation inc»euses with harmonic number 
and strength parameter. In addition, Herma i Seow rediet ; 
except qualitatively, the phase relations in the distorted 


wave. 
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APPENDIX A 


Tite Oe kOGkAM PARAMETERS AND OPTIONS 


The purpose of this program was to convert up to four 
channels of analog data into a series of digital points on 
seven-track magnetic tape. The program had several input 
parameters which had to be provided to the computer via 
Seemee corinter arter the PORTRAN program had been read into 
the system, 

NSAMP (input 1.) was an integer which corresponded to the 
number of data or samples to be taken from each channel. 

For this investigation a sampling rate of 10 kHz was used 
and NSAMP was chosen to span one cycle of the input waveform. 

NCHAN (input 2.) was an integer which specified the 
number of channels of analog data being converted. 

MeEee (input 2.) specified — many records were to be 
written on tape each time sampling was initiated. A record 
was defined as the data set corresponding to the product of 
NCHAN and NSAMP . BeeaMcesOrechnemhigh Sampling rate used ie 
was necessary to set this value to l. 

ITAPE (input 4.) designated which magnetic-tape trans- 
port held the seven-track tape. The integer used here had 
EOwec@nmrespond Ehat selected on the “unit select" switch on 
the particular tape unit used. 

NDEL (input 5.) determined how rapidly the data on the 


Seven-track tape was read back when using program option 


a2 





An example of the input parameters as typed on the 


teleprinter might be 


Veta ont, 1000, NCHAN = 3, ITAPE = 1, 


NDEL 


20* (carriage return). 


Mmiewconterol! symbol * indicated to the computer that the 
input parameters had been set. Care was taken when choos- 
mageNoAMP and NCHAN to insure that their product did not 
exceed the first dimension of IBUF found in the first FORTRAN 
Program statement. 

When the parameter list had been given the teleprinter 


responded with 
"OPTION = (I1)". 


MeeweOmpulcer was then given a one-digit number corresponding 
Pemene list OL available options. 

RESTART (option 1.) was used whenever any one or all of 
the input parameters was to be changed. 

GO (option 2.) was used when analog to digital conver- 
sion was to be done. Upon receiving this option the compu- 
ter awaited the initiating trigger from switch DS-1l to 
start the actual conversion process. Seria (renege freon, Sa - 
caused one record of data to be written on tape. When a 
sufficient number of records were written, DS-l was 
Penge reamwith DS-2 in the “up" position. This caused the 


last record to be written followed by a teleprinter response 


Orr LON = (11)”. 
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Hi sOntee tie Optlon 3.) was used to put a mark on the 
tape to separate or isolate groups of records. Groups of 
records separated by END OF FILE marks are called files. 

SKIPFILES (option 5.) were used when files already on 
tape were to be skipped over. When used, a four-digit 
number was given the computer to specify the number of 
files to be skipped. 

LIST ONE RECORD (option 6.) was used in reading data 
from the seven-track tape. Use of this option caused data 
Piemeren On tape to be written on the page printer after 
faeces peciiied the number of data to print. 

Emor ONE RECORD (option 7.), when called, Saeed one 
record to be read from the tape, converted to an analog 
Bee reeanecmread) 2O terminals T420 through F422 cn the analog 
board feOrmesponding tO Channels 1 through 4). The signal 
could then be monitored on various recording devices. 
Because of the short duration of the Signal it was best 
displayed on the available pen recorder whose input terminals 
are found on the analog board (Figure 3). The time base 
of the signal returned from the tape could be lengthened or 


shortened by varying the value of the input parameter NDEL. 
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